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Abstract 

We report first measurements of the doubly charmed baryonic B decays B — > A+A~K. The 
B + — > A+A~K + decay is observed with a branching fraction of (6.5lJ;g ± 1.1 ± 3.4) x 1CT 4 and a 
statistical significance of 15.4 a. The B° — ► A+A~K° decay is observed with a branching fraction 
of (7.9I2 3 ^ =t 4.1) x 10~ 4 and a statistical significance of 6.6 a. The branching fraction errors 
are statistical, systematic, and the error resulting from the uncertainty of the — ► pK~ir + decay 
branching fraction. The analysis is based on 357 fb -1 of data accumulated at the T(45) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e _ collider. 

PACS numbers: 14.20.Lq, 14.40.Nd 
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Recent 
reported 



ajiumber of studies of single charmed baryon production in B decays have been 
The measured branching fractions of the two-body single charmed 




baryon decays B° — > Afv La] and B — > S°(2455)p [4] are significantly smaller than theoret- 
ical expectations p, la, ullsj]. The multi-body single charmed baryon decays B — > A+p7r(7r) 
were found to have branching fractions about one order of magnitude larger than the cor- 
responding two-body decays, but still below theoretical predictions. While single charm 
production proceeds via a b—>cud quark transition, production of two charmed particles 
occurs via a 5— »ccs transition. In contrast to the single charmed baryon production the 



two-body doubly charmed baryon B decay B + — > S°A+ £ | recently observed at Belle has a 
branching fraction comparable to theoretical predictions 5J]. It would be interesting to check 
whether theory can describe multibody double charmed decays. In this paper we report the 
first observation of the B + — ► A+A~K + and B° — > A+A~K° decays, which are three-body 
decays that proceed via a 5— >ccs transition. Inclusion of charge conjugate states is implicit 
unless otherwise stated. The analysis is based on a data sample of 357 fb -1 accumulated 
at the T(45) resonance with the Belle detector at the KEKB asymmetric-energy collider 
corresponding to 386 x 10 6 BB pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K® mesons and to identify muons 
(KLM). The Belle detector is described in detail elsewhere |lOj. Two different inner detector 
configurations were used. For the first sample of 152 x 10 6 BB pairs, a 2.0 cm radius beam- 
pipe and a 3-layer silicon vertex detector were used; for the latter 234 x 10 6 BB pairs, a 
1.5 cm radius beam-pipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
used (llj. We use a GEANT based Monte Carlo (MC) simulation to model the response of 
the detector and determine its acceptance jl^ . 

We detect the A+ via the A+ — > pK~TT + , pK° and Att + decay channels. When a A+ and 
A~ are combined as B decay daughters, at least one of A^ is required to have been recon- 
structed via the pK^ir^ decay process. For each charged track, the particle identification 
(PID) information from the CDC, ACC and TOF is used to construct likelihood functions 
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Lp, Lx and L n for the proton, kaon and pion assignments, respectively. Likelihood ratios 
L a /(L a + L b ) are required to be greater than 0.6 to identify a particle as type a, where b 
denotes the other two possible hadron assignments from the three possiblities: proton, kaon 
and pion. For the main mode B + — > A^A~K + , A+ — > pK~ir + , A~ — > pK + ix~ the PID effi- 
ciency for the primary K + is about 95%. Efficiencies for protons, kaons and pions from A+ 
decays are about 98%. The misidentification probability for pions (or kaons) to be identified 
as kaons (or pions) is less than 5%. The probability for pions or kaons to be identified as 
protons is less than 2%. Tracks consistent with an electron or muon hypothesis are rejected. 
A A+ candidate is selected if the mass of its decay products is within 0.010 GeV/c 2 (2.5 a) 
of the nominal A+ mass. 

Neutral kaons are reconstructed in the Kg— >tt + 7t~ decay. Candidate A baryons are re- 
constructed in the decay A— >p7r~. We apply vertex and mass constrained fits for the K° 
and A candidates to improve the momentum resolution. The intersection point of the K° 
and A candidate daughter tracks must be displaced from the beam interaction point (IP): 
the flight distance should be more than 0.5 mm. A K° candidate is selected if the mass of 
its decay products is within 7.5MeV/c 2 (3cr) of the K° mass. A A candidate is selected if 
the mass of its decay products is within 2.5 MeV/c 2 (2.5 a) of the A mass. 

The B candidates are identified using the beam-energy constrained mass My, c and 
the mass difference AMg. The beam-energy constrained mass is defined as Mt> c = 
\/-^bcam ~~ (J2Pi) 2 , where -Ebeam is the beam energy, and pi are the three-momenta of the B 
meson decay products, all defined in the center-of-mass system (CMS) of the e + e~ collision. 
The mass difference is defined as AMb = M(B) -mg, where M(B) is the reconstructed 
mass of the B candidate and me is the world average B meson mass. The parameter AMb 
is used instead of the energy difference AE = (52 Ei) ~ E beiim , where E^ is the CMS en ergy 
of the B decay products, since AE shows a correlation with M^ c , while AMb does not jl3j |. 
M(B) = ^/E(By-(J2p t ) 2 , where E(B) = E(A+) + E(A-) + E(K), E(A+) = ^ A 2 ++m* c+ , 
p A + is the A+ momentum measured via its decay products and m A + is the value of the A+ 
baryon mass Q- We select events with M bc > 5.20 GeV/c 2 and \AM B \ < 0.20 GeV/c 2 . The 
prompt K + or the reconstructed Kg trajectory and the A+/A~ trajectories are required to 
form a common B decay vertex. If there are multiple candidates in an event, the candidate 
with the best Xb f° r the B vertex fit is selected. The B vertex fit is performed without 
additional mass constraints for known particles. 
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FIG. 1: Candidate (a, b) B+ 



A+K K+ and (c, d) B° 

distribution for M hc > 5.27GeV/c 2 and (b, d) M bc distribution for |AM B | < 0.015 GeV/c 2 . Curves 
indicate the fit results. 



A+A c K° decay events: (a, c) AM B 
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Figure ^ shows AMb and M^ c projections for selected B + — > A+Ajii" + and 
5° — ► A+A~i^° decay events. The AMg projection is shown for M^ c > 5.27GeV/c 2 and 
the Mbc projection is shown for |AM#| < 0.015 GeV/c 2 . The widths determined from single 
Gaussian fits to MC generated events are 2.7 MeV/c 2 and 3.3 MeV/c 2 for M^ c and AMb, 
respectively. A two-dimensional binned maximum likelihood fit is performed to determine 
the signal yield. The AMb distribution is approximated by a Gaussian for the signal plus 
a first order polynomial for the background, and the Mh r distribution is represented by a 



single Gaussian for the signal plus an ARGUS function [15j for the background. The signal 
shape parameters are fixed to the values obtained from a fit to a MC simulation. All yields 
and background shape parameters are allowed to float. 

From the fit we obtain signal yields of 48.5l 6 ;g and 10.5i|;f events with statistical sig- 
nificances of 15.4cr and 6.6a, for B + — > A+A~K + and B° — > A+A~K°, respectively. The 
significance is calculated as 



-2 ln(£ /£ max ), where £ max and C denote the maximum 
likelihoods with the fitted signal yield and with the yield fixed at zero, respectively. 

The branching fraction Bij for the i-th A+ decay and the j-th A~ decay mode are cal- 
culated as = Nij /[N BB SijBi(A+)Bj(A~)], where is the B signal yield. The de- 
tection efficiencies, determined from MC simulation. The A+ decay branching 
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fractions £>j(A + ) are converted to the product £>(A+ — > pK 7i + )Ti/T(pK n + ) to isolate 
the common uncertainty from the branching fraction of A+ — > pK~n + . The values of 
Ti/T(pK-7r + ) are (0.47 ± 0.04) and (0.180 ± 0.032) for pK° and Att+ modes, respec- 
tively Q|. The overall detection efficiency e for the total signal yield N is calculated 
as J2 e ij[Fi/T(pK~TT + )][Tj/r(pK~'n- + )]. The overall branching fraction is calculated as 
N s /[N B beB(A+ — > pK~n + ) 2 ], using the overall signal yield As and the decay branching 
fraction £>(A+ — > pK~ir + ) = (5.0 ± 1.3)% The detection efficiencies are calculated to 
be 7.79% for the B+ -> A+A~if + decay and 1.38% for the 5° -> A+A~K° deca y- 

The number of pairs A" B £ is (386 ±4) x 10 6 . The fractions of charged and neutral 
B mesons are assumed to be equal. We obtain branching fractions of 



B{B + -> A+A-^ 



(6.5±J;g ±1.1 ±3.4) x 10" 4 and 



B(B° -> A+A-iT ) = (7.9±|;| ± 1.2 ± 4.1) x 10" 



where the first and the second errors are statistical and systematic, respectively. The last 
error is due to the 52% uncertainty in the absolute branching fraction, £>(A+ — > pK~7i + ). 

Systematic uncertainties in the detection efficiencies arise from the track reconstruction 
efficiency (8% - 10% depending on the process, assuming a correlated systematic error of 
about 1% per charged track); the PID efficiency (9% - 10% assuming a correlated systematic 
error of 2% per proton and 1% per pion or kaon); 3-body decay model uncertainty (11% for 
the B + -> A+A~K+ decay and 5% for the B° -> A+A^fT decay); and MC statistics (1% 
- 2%). The other uncertainties are associated with T(A^)/T(pK^n + ) (2% - 3%); and the 
number of N B g events (1%). The total systematic error is 17% for B + — > A+A~K + and 
15% for B+ -> A+A-K . 
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FIG. 2: M(A±) mass distributions for (a) B + -> A+A C 7K+ and (b) 5° -> A+A^K decay candi- 
dates in the i? signal region. Curves indicate the fit results. 
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Figure El shows the mass distributions M(Af) for B candidates in the signal region 
\AM B \ < 0.015 GeV/c 2 and M bc > 5.27GeV/c 2 . The M(A±) mass distributions are shown 
for \M(AJ) — m A +\ < 0.010 GeV/c 2 . The curves show the results of a fit with the sum of 
a Gaussian and a linear background. The means and widths of the Gaussians are fixed to 
values obtained from fits to MC samples. For B + — > A+ A~K + decay, we obtain a A+ yield of 
39.5l£jj events and a A~ yield of 48.2±£ ' events. For B° -> A+A~K°, yields of 11.4±|;| and 
10.0 jzllf events are obtained from the A+ and A~ distributions, respectively. These values 
are consistent with the B signal yields given above. 

We consider possible contributions from other B decays, which could give a B signal in 
the AE and AMb distributions, but should produce a uniform distribution in the A+ mass 
region. To assess this type of background, we analyze the A+ sideband 0.015 GeV/c 2 < 
|M(A+) - m A+ | < 0.055 GeV/c 2 and |M(A~) - m A -\ < 0.010 GeV/c 2 , and A" sideband 



0.015 GeV/c 2 < |M(A" 



< 0.055 GeV/c 2 and |M(A+) - m A+ \ < 0.010 GeV/c 2 . 



We conclude that other B decays contribute less than 1.7 events at 90% C.L. in the 



B- 



A+A~K + mode and less than 0.2 events at 90% C.L. in B° -> AfA~K°; both con- 



tributions are neglected. 
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FIG. 3: M(A+A~) mass distributions for (a) B A 
candidates in the B signal region. Points with error bars are data. Open histograms - MC for a 
uniform phase space distribution. Hatched histograms - normalized A+ sideband data. 



Figure El shows the M(A+A~) mass distributions for (a) B + — > A+A~K + decay candidates 
and (b) B° — > A+A~i"f° decay candidates in the B signal region |AMb| < 0.015 GeV/c 2 and 
M bc > 5.27GeV/c 2 . No deviations from phase space distributions are evident. 

In summary, we have reported the first measurement of the doubly charmed baryonic 
B decay B + — > A+A~K + with a branching fraction of (6.51q 9 ± 11 ± 3.4) x 10~ 4 and a 
statistical significance of 15.4 a, and the B° — > A+A^i^ decay with a branching fraction of 
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(7-9-2 3 ^ 1-2 i 4.1) x 10~ 4 and a statistical significance of 6.6 a. These three-body doubly 
charmed B decay branching fractions are about the same order of magnitude (or slightly 
smaller) than the branching fraction of the two-body doubly charmed decay B + — > S°A+, 
which is due to the same b— >ccs quark transition, also observed by Belle [jj. The behavior of 
these b—>ccs decays is qualitatively different from single charmed baryon decays, where three- 
body decays have bigger branching fractions than two-body decays. The obtained branching 
fraction is by five to six orders of magnitude higher than expected from naive estimation 
for the B — > A^A~K decay with color suppression, which is also highly suppressed by phase 
space All this needs further experimental and theoretical study. 
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